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Abstract 

Nuclear shadowing in DIS at moderately small x is suppressed by the nuclear 
formfactor and depends on the effective mass of a hadronic fluctuation of the virtual 
photon. We propose a solution to the problems (i) of how to combine a definite 
transverse size of the fluctuation with a definite effective mass, and (ii) of how to 
include the nuclear formfactor in the higher multiple scattering terms. Comparison of 
the numerical results with known approximations shows a substantial difference. 



1. Introduction 

Shadowing in deep-inelastic scattering (DIS) off nuclei is a hot topic for the last two 
decades. In the infinite momentum frame of the nucleus it can be interpreted as a result 
of parton fusion leading to a diminishing parton density at low Bjorken x [Q - {§]. A more 
intuitive picture arises in the rest frame of the nucleus where the same phenomenon looks like 
nuclear shadowing of hadronic fluctuations of the virtual photon || - [|T!|. To crystallize 
the problem and its solution we restrict ourselves in this paper to only quark-antiquark 
fluctuations of the photon, neglecting those higher Fock components which contain gluons 
and qq pairs from the sea. The lifetime of the qq fluctuation (called coherence time) is given 
by 

tc ~ Q 2 + M 2 { ) 

where v is the photon energy, Q 2 its virtuality and M is the effective mass of the qq pair. 
Provided that the coherence time is much longer than the nuclear radius, l c ^> Ra, the 



total cross section on a nucleus reads 13 



°L A (x,Q 2 ) = 2 fd 2 bj d 2 rG r (Q 2 ,r) jl - exp [_I a (r)T(&)] } 

= 2 Jd 2 b {l-(exp[-ia(r)T(6)])} . (2) 

Here G 7 * (Q 2 , r) characterizes the probability for the photon to develop a qq fluctuation with 
transverse separation r. The condition t c 3> Ra insures that the r does not vary during 
propagation through the nucleus (Lorentz time dilation). Then the qq pair with a definite 
transverse separation is an eigenstate of the interaction with the eigenvalue of the total cross 
section a{r). Therefore, one can apply the eikonal expression (0) for the interaction with 
the nucleus. The nuclear thickness function T(b) = j*^ dz pA(b, z) is the integral of nuclear 
density over longitudinal coordinate z and depends on the impact parameter b. 

The color dipole cross section a(r) introduced in |13[ vanishes like r 2 at small r — > due 
to color screening. This is the heart of the phenomenon called nowadays color transparency 
|14] , [L3|, [15]. For this reason nuclear shadowing in (§) is dominated by large size fluctuations 
corresponding to highly asymmetric sharing of the longitudinal momentum carried by the 



q and q |16j, [F| ^ ^2| . This leads to Q 2 scaling of shadowing. 

Note that the averaging of the whole exponential in makes this expression different 
from the Glauber eikonal approximation where er(r) is averaged in the exponent. The 
difference is known as Gribov's inelastic corrections |T7|]. In the case of DIS the Glauber 
approximation does not make sense, and the whole cross section is due to the inelastic 
shadowing. 

For the other case, t c ~ Ra, one has to take into account the variation of r during the 
propagation of the qq fluctuation through the nucleus. At present this can only be done for 



the double scattering term |L2] in the expansion of the exponential in (g), 

(T) Jd 2 bF 2 (q,b) + (3) 



*& N ~ 4 (a{r)) 
or in hadronic representation |L8] , 



1- WfrtfdM 



where the mean nuclear thickness and the formfactor read, 



&)+..., (4) 

4=0 



(T) = j J d%T\b) , (5) 
1 °° 

F A (q,b) = -- J dzp A (b,z)e^ , (6) 

— oo 

with longitudinal momentum transfer g = l/t c given by (jl]). In the case of (|3]) the uncertain 
fluctuation mass is fixed at M 2 = Q 2 , q = 2m^x. Two expressions (||) and (f|) are related 
since the integrated forward diffractive dissociation cross section 7*iV — > XiV equals to 
(a 2 )/167T. 

There are two problems remaining which are under discussion: 

• How the nuclear formfactor can be included in the higher order scattering terms which 
are of great importance for heavy nuclei? For instance, the shadowing term in (|3]), 
(|]) for lead is of the order of one at low x, so the need of the higher order terms is 
obvious. 
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• Even for the double scattering term in (Q) it is still unclear which argument should 
enter the formfactor. Indeed, the effective mass of the qq fluctuation needed for the 
coherence time in (|1|) cannot be defined in the quark representation with a definite qq 
separation. On the other hand, Eq. ([|) exhibits an explicit dependence on Mx and 
the longitudinal momentum transfer is known. However, unknown in this case is the 
absorptive cross section of the intermediate state X. 

We suggest a solution of both problems in the next section. The goal of this paper is 
restricted to the study of the difference between the predictions of the correct quantum- 
mechanical treatment of nuclear shadowing and known approximations. We do it on an 
example of the valence qq part of the photon and neglect the higher Fock components 
containing gluons and sea quarks, which may be important if to compare with data especially 
at very low x. Nuclear ant i- shadowing effect is omitted as well, since we believe it is beyond 
the shadowing dynamics (e.g. bound nucleon swelling). Numerical results and a comparison 
with the standard approach are presented in section 3. 

2. The Green function of a qq pair in nuclear medium 

We start with the generalizing of eq. (fj) for the case l c < Ra, 

i 

olot A (x,Q 2 ) = J " d 2 b J daa? ot A (XiQ 2 ;b,a) , (7) 
o 

where 

a£ A (x,Q 2 ;b,a) = T(b) J d 2 r \%*(r, a)\ 2 a{r) 

oo oo 

- 2Rey dz x p A (b,z x ) J dz 2 p A (b, z 2 ) A(z u z 2 , a) . (8) 

— OO 21 

The first term in r.h.s. of @ corresponds to the second, lowest order in a(r)T(b), term 
in expansion of the exponential in @. The shadowing terms are contained in the second 
term in (H). (r, a) is the (non- normalized) wave function of the qq fluctuation of the 
virtual photon, where a is the fraction of the light-cone momentum of the photon carried 
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by the quark. An explicit expression of transverse and longitudinally polarized photons can 
be found in 

The function A(zi,Z2,oc) in (|j) reads, 
1 f 

A{z u z 2 ,a) = - / d 2 r 1 d 2 r 2 ^(f 2 ,a)W(f 2l z 2 ;r 1: z 1 )^ r (f 1 ,a)a(r 2 )a(r 1 )e iq ^ {z2 ~ Zl) , 



with 



Q 2 a(l — a) +m 2 q 
2va(l — a) 



(9) 
(10) 



This expression was first suggested in unpublished paper p0 



The second (shadowing) term in (§D is illustrated in fig. 1. At the point Z\ the photon 




Figure 1: A cartoon for the shadowing (negative) term in ([|). The Green 
function W(r 2 , z 2 ;ri,zi) results from the summation over different paths 
of the qq pair propagation through the nucleus. 



diffractively produces the qq pair (7*iV — > qqN) with transverse separation ri. The pair 
propagates through the nucleus along arbitrarily curved trajectories (should be summed 
over) and arrives at the point z 2 with a separation f 2 . The initial and the final separations 
are controlled by the distribution amplitude \E f 7 *(r). While passing the nucleus the qq 
pair interacts with bound nucleons via the cross section <j(r) which depends on the local 
separation r. The function W(f 2 , z 2 ;fi,Zi) describing the propagation of the pair from z\ to 

5 



z 2 also includes that part of the phase shift between the initial and the final photons, which 
is due to transverse motion of the quarks, while longitudinal motion is already included in 
( |T0| ) via the exponential. 

Thus, Eq. (^|) does not suffer from either of the two problems of the approximations (0) 
- (p. The longitudinal momentum transfer is known and all the multiple interactions are 
included. 

The propagation function W(f 2 , z 2 \ fx, z\) in @ satisfies the equation p0| , 

. dW(f 2 ,z 2 ;f 1 ,z 1 ) A(r 2 ) 



dz 2 



2ua(l — a 



W(f 2 ,z 2 ;f 1 ,z 1 ) 



vfa) PaO>, z 2 ) W(f 2 , z 2 ; n, z\) 



with the boundary condition W(f 2 , z±; fi, z\) = 5(f 2 — fi). The Laplacian A(r 2 ) acts on 
the coordinate f 2 . The full derivation of ([11]) will be given elsewhere. Here we only notice 
that it looks natural like Schrodinger equation with the kinetic term A/[2i/a(l — a)] which 
takes care of the varying effective mass of the qq pair and provides a proper phase shift, 
and z 2 plays the role of the time. The imaginary part of the optical potential describes the 
absorptive process. 

In the "frozen" limit v — ► oo the kinetic term in flTT| ) can be neglected and 



W(f 2 , z 2 ; f u zx) = S(r 2 - fi) exp 



1 

-cr(r 2 ) J dzp A (b,z) 

21 



and with q r , 



(12) 

one arrives at result 



When this expression is substituted into (||) - 
d) with G r {Q 2 ,r) = f^da |^ 7 *(f,a)| 2 . 

We can also recover the approximation (|3[) - (f|) if one neglects the absorption of the qq 
pair in the medium. Then W becomes the Green function of a free motion, 



W{r 2 ,z 2 \r x ,z x )\ a _ 



2tt 



d 2 k exp 



ik(f 2 — fi) + 



ik 2 (z 2 - zi) 
2va{\ — a) 



(13) 



where k is the transverse momentum of the quark. 

With this expression the shadowing term in (|$D reproduces the second term in (f|). 
Indeed, the amplitude of the photon diffractive dissociation in the plane wave approximation 



reads, 

f dd (k) = 

Therefore, ® can be represented as 



^ f d 2 r (r, a) cr(r) e 



'kr 



A(z 1 ,z 2 ,a) = ^- I d 2 k \fdd(k)\ z exp 



Q 2 a{l - a) + m 2 + k 21 
(z 2 -z{)- 



(14) 



(15) 



2va{\ — a) 

Taking into account that M\ = {m 2 + k 2 )/a{l — a) is the effective mass squared of the qq 
pair and substituting (|i~5|) to (§) we arrive at eq. (|). 



3. Numerical results 

We calculate nuclear shadowing for calcium and lead from the above displayed equations. 
As was mentioned in the Introduction, only the valence gg-part of the photon is taken into 
account, but the higher Fock components containing gluons and sea quarks are neglected, 
as well as the effect of anti-shadowing. Therefore, we do not compare our results with data, 
but only to the standard approach (|3|) - (f|). 



We do the same calculations again, using the free Green function ([13]). This makes it 
possible to disentangle between the influence of higher scattering terms and the formfactor. 

We approximate the cross section by the dipole form <r(r) = Cr 2 , C ~ 3, which is a good 
approximation at r > 0.2—0.3 fm ||21|| . However, we calculated the proton structure function 
F 2 (x,Q 2 ) perturbatively (we fixed the quark masses at m q = 0.3 GeV, m s = 0.45 GeV and 
m c = 1.5 GeV) what leads to an additional logarithmic r-dependence at small r. This is 
important since results in the double- log Q 2 dependence of F 2 . Nuclear shadowing, however, 



is dominated by soft fluctuations with large separation |fL2| , therefore, the dipole form of 
the cross section is sufficiently accurate. 

We use a uniform density for all nuclei, pa = 0.16 fm~ 3 , what is sufficient for our 
purpose, comparison with the standard approach calculated under the same assumption. 

Within these approximations it is possible to solve (|TTj) analytically. The solution is the 
harmonic oscillator Green function with a complex frequency pT| , 

2f 2 ■ fi 



W (f2,z 2 ;r u z 1 ) 



27rsinh(cuA^ GXP I 2 



( r 2 + r i) c °th (^A; 



sinh (uAz) 



, (16) 



where 



Az = z 2 — Z\ 

2 ■ C PA 



va (1 — a) 

a 2 = —i Cpavol (1 — a) . (17) 

This formal solution properly accounts for all multiple scatterings and finite lifetime of 
hadronic fluctuations of the photon, as well as for fluctuations of the transverse separation 
of the qq pair. 

The results of calculations are shown in fig. The dashed curves show predictions of 
(0) which we call standard approach. The mean values of a 2 and a are calculated using the 
same qq distribution functions Jl9], |J as in (H) and the intermediate state mass is fixed at 
M 2 = Q 2 . At low x < 0.01 shadowing saturates because q = Im^x <C 1/Ra- The thin solid 
curve also corresponds to a double scattering approximation, i.e. absorption (the second 
term in (|TT|) ) is omitted. However, the formfactor is treated properly, i.e. the kinetic term in 
( jljD taking into account the relative transverse motion of the qq pair, correctly reproduces 
the phase shift. The difference between the curves is substantial. The thin solid curve does 
not show saturation even at x — 0.001. 

The next step is to do the full calculations and study importance of the higher order 



rescattering terms in (11). The results are shown by the thick solid curves. Higher or- 
der scattering brings another substantial deviation (especially for lead) from the standard 
approach. At very low x the curves saturate at the level given by @). 

4. Conclusions and outlook 

We suggest a solution for the problem of nuclear shadowing in DIS with correct quantum- 
mechanical treatment of multiple interaction of the virtual photon fluctuations and of the 
nuclear formfactor. We perform numerical calculations for qq fluctuations of the photon 
and find a significant difference with known approximations. Realistic calculations to be 
compared with data on nuclear shadowing should incorporate the higher Fock components 
which include gluons. The same path integral technique can be applied in this case. The 
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x- dependence of the dipole cross section a(r, x) (correlated with r p2| ) should be taken into 
account. One should also include the effect of ant i- shadowing, although it is only a few 
percent. A realistic form for the nuclear density should be used (this can be done replacing 



F 8 , A /AF 8 , N 




Figure 2: Nuclear shadowing for calcium and lead. The dotted curve is calculated in 
the standard approach ([|). The thin solid curve corresponds to the double scattering 
approximation with the free Green function, and the thick solid curve shows the full 
calculation, ftfflj). 
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PA{b,z) by a multistep function like in We are going to settle these problems in a 

forthcoming paper. 
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